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CP-odd phase of the stop trilinear soft breaking term, aig{At), are studied. We discuss the CP- 
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^ " CP-nonconserving effects in the reaction pp ^ tb + X ^ W~^bb + X, driven by the supersymmetric 
Oh 

o I nonconserving effects in both production and the associated decay amplitudes of the top. We find 
^ that, within a plausible low energy scenario of the MSSM and keeping the neutron electric dipole 
j> [ moment below its current limit, a CP-violating cross-section asymmetry as large as 2 — 3% can 
^ . arise if some of the parameters lie in a favorable range. A partial rate asymmetry originating only 
• in the top decay t — > W'^b is found to be, in general, below the 0.1% level which is somewhat 
smaller than previous claims. For a low tan/3 of order one the decay asymmetry can reach at the 
most ~ 0.3%. This (few) percent level overall CP-violating signal in pp ^ tb + X ^ W'^bb + X 
might be within the reach of the future 2(4) TeV pp Tevatron collider that may be able to produce 
~ 10000(~ 30000) such tb events with an integrated luminosity of 30 fb~^. In particular, it may 
be used to place an upper bound on arg(Aj) if indeed arg(yu) 0, as implied from the present 
experimental limit on the neutron electric dipole moment. The partial rate asymmetry in the top 
decay ( ~ few x 10^^ ) may also be within the reach of the LHC with ~ 10'' pairs of tt produced, 
provided detector systematics are sufficiently small. We also show that if the GUT-scale universality 
of the soft breaking trilinear A terms is relaxed, then the phases associated with arg Au and arg Ad 
can take values up to ~ few x 10~^ even with squarks and gluino masses of several hundred GeV's 
without contradicting the experimental limit on the neutron electric dipole moment. 



^Address after August 1, 1997. 
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1. Introduction 



In the Tevatron pp collider, top quarks will be mainly produced as pairs of ti via an s-channel 
gluon exchange. Nonetheless, the subleading Electro Weak (EW) production mechanism of a single 
top, forms a significant fraction of the tt pair production. It will therefore be closely scrutinized in 
the next runs of the Tevatron In particular, the production rate of tb (and the charged conjugate 
pair) through an s-channel ly-boson, pp —>■ W ^ tb + X , (the corresponding partonic reaction is 
ud ^ W tb) is about 10% of the ti production rate Throughout this paper we will always 
refer only to the s-channel W exchange when discussing the reaction pp ^ tb + X. 

In a previous work, 0, we have shown that a new CP- violating phase from a Two Higgs Doublet 
Model (2HDM) of type II can give rise to CP-violating asymmetries of the order of a few percent in 
the reaction pp —* tb + X. In that work we had considered one example of a supersymmetric (SUSY) 
mediated, one-loop, CP- violating diagram (gluino exchange diagram). The one SUSY diagram that 
was considered in [H was found to be p-wave suppressed near threshold, thus yielding CP asymmetries 
of the order of 0.1%. Clearly a complete study of the possible SUSY-CP effect is required in order 
to be able to estimate the true expected magnitude of a potential CP-nonconserving effect that is 
driven by new CP-odd phases of low energy SUSY dynamics. 

Therefore, in this paper we wish to extend our previous study and explore the complete CP- 
violating effect in the process pp tb + X W~^bb + X, induced specifically by the complex soft 
trilinear parameter associated with the superpartner of the top. At, in the Minimal Supersymmetric 
Standard Model (MSSM). We will separately discuss the SUSY CP-nonconserving effect in produc- 
tion and decay amplitudes of the top and show how these two asymmetries combine to give the 
overall asymmetry in pp ^ tb + X ^ W~^bb + X. 

The experimental limit on the Neutron Electric Dipole Moment (NEDM), (i„ < 1.1 x 10^^^ e-cm 
[0, places a severe constraint on the phase of the SUSY Higgs mass parameter /i. In particular, 
arg(/i) < O(10~^) m, H, H, 01 is essentially inevitable for a typical SUSY mass scale of the order 
of a few hundreds GeV. Therefore, the only significant SUSY-CP-odd phase, that can potentially 
drive notable CP-nonconserving effects in top quark reactions, is the phase in the stop soft trilinear 
breaking term At. As will be shown below, CP- violating effects that can arise from the other Ag- 
terms (associated with the lighter quarks g^) are extremely suppressed being proportional to the 
masses of the light quarks. In particular, when m^^ — there is no mixing between the left and 
the right squarks. The two mass eigenstates of the supersymmetric partners of the light quarks are 
therefore expected to be nearly degenerate, thus not playing any role in CP-violation effects at high 
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enough energies. Of course, this is not the case for the NEDM which is driven by the shght deviation 
from degeneracy of the supersymmetric partners of the u and the d quarks when arg(/i) = 0. 

Moreover, if the phases of A^, and At are not correlated at the EW scale, which is indeed 
possible if the GUT-scale universality of the A terms is relaxed, then the experimental limit on 
the NEDM cannot put any further constraint on a.Tg{At). Therefore, it is extremely important 
to explore other avenues for constraining a,Tg{At). Of course, the most natural place to look for 
CP-nonconserving effects, which are driven by arg(Af), is high energy reactions involving the top 
quark. Thus, in the limit arg(yu) — > 0, the only important phenomenological CP-violating SUSY 
parameter, that enters any CP-odd effect at high energies, resides in the imaginary part of the 
ti — tn mixing matrix elements lm{Zl^* Z^"^) (see appendix A). Being proportional to a,Tg{At), it 
vanishes for a.Tg{At) 0. Besides, aig{At) may play an important role in explaining the observed 
baryon asymmetry in the universe. Recently, it was claimed that with arg(/i) — 0, t squarks can 
mediate the charge transport mechanism needed to generate the observed baryon asymmetry even 
with squark masses ~ few x GeV, provided that a.Tg{At) is not much suppressed 

In this paper we discuss the CP-odd effects, induced by lm{Z}" Z^^) oc arg(Ai) in the reaction 
pp ^ tb + X and in the subsequent top decay t — > W~^b. We find that, with maximum CP-violation 
(i.e., \lm{Z}'^* Zt^)\ = 1/2), a CP-violating asymmetry in the cross-section can reach 3% for squark 
masses at around 0.5 TeV and a light stop mass below 100 GeV. We also find that, in the same range 
of squark masses, a Partial Rate Asymmetry (PRA) effect int—* W^h is predominantly below 0.1% 
and can be slightly above 0.1% for low tan/5. In a narrow window of the SUSY parameter space, 
with tan (3 ~ 0{1), it can reach at the most ~ 0.3%. This is somewhat less than the estimates in the 
literature of a few percent asymmetry in top decay in the MSSM 0. Therefore, with a predicted few 
percent overall asymmetry in the cross-section and decay, the reaction pp ^ tb + X ^ W^hb + X 
can serve to constrain aig{At) in the future runs of the upgraded Tevatron at Fermilab. 

Also, as a by product of our analysis, we found that when the Higgs mass term, /x, is real at 
the EW-scale (or has a very small phase) and bearing theoretical uncertainties associated with the 
naive quark model approach (see e.g., 0]), the experimental limit on the NEDM can be naturally 
accommodated with phases of the soft trilinear breaking A terms at the order of ~ few x 10~^ even 
with squark and gluino masses below 500 GeV. Therefore, we believe that there is no compelling 
argument for the SUSY CP-odd phases to be necessarily less than or of the order of 10~^ for squarks 
masses below 1 TeV, contrary to previously and commonly used claims. 

We wish to emphasize that the analytical formulations given in this paper hold for the most 
general low energy realization of the MSSM. Moreover, the PRA effect in the top decay t —>■ W~^b, 
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presented in this paper, does not depend on the specific production mechanism of the t and t. 
Therefore, the LHC which may be able to produce ~ 10^ pairs of tt may, in principle, be able to 
detect a PRA effect in t —>■ W^h even as small as ~ few x 10"'^, provided detector systematics are 
sufficiently small. 

The paper is organized as follows: in section 2 we construct a plausible phenomenological low 
energy SUSY spectrum. In section 3 we present a complete, self contained analytical derivation of 
the CP asymmetry in both production and decay processes of the top in a general MSSM framework, 
keeping all possible SUSY CP-odd phases. In section 4 we discuss our choice of the low energy CP 
phases and its implications on the NEDM. In section 5 we present our numerical results in the 
limit arg(/i) = and in section 6 we summarize. The relevant SUSY lagrangian pieces along with a 
detailed derivation of the various SUSY mixing matrix elements is given in appendix A. In appendix 
B we list the SUSY CP-violating phases, in the general MSSM scenario, that are responsible for the 
CP-asymmetry in our reaction while a description of the one-loop integrals is given in appendix C. 

2. Low Energy MSSM Phenomenology 

The most general low energy = 1 minimal supergravity (SUGRA) SU{'i) x SU(2) x U{1) 
invariant lagrangian (which apart from gravitational interactions, is essentially identical at low en- 
ergies to a theory with softly broken supersymmetry), that consists of three generations of quarks, 
two Higgs doublets and the SU{3) x SU{2) x f/(l) gauge fields, along with their supersymmetric 
partners, can be written as [|1T], |I^: 



eij is the antisymmetric tensor with eu = 1 and the usual convention was used for the superfields 




(1) 



where W is the superpotential and is given by: 



(2) 



Q, U, L, R and H |jTT|. The capital index indicates the generation (i.e., J, J = 1, 2 or 3). 



£soft consists of the soft breaking terms and can be divided into three pieces: 



•gaugino 




(3) 
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These are the soft-supersymmetry breaking gaugino and scalar mass terms and the trihnear couphng 
breaking terms. In particular they are given by [p!3| : 



^gaugino = ^{mi\B\B + "^2A^A^ + msAcAc) , (4) 

^scalar = -m^J i// | ' - m^J ff^ 1 2 _ ^2 | 1 2 _ 

ml\R\' - mllQf - ml\D\' - ml\U\' , (5) 



Arilinear = eij{gu AuQ'' + QdAdQ^ HID + 

gEAEL'HiR + iiBH\Hi) . (6) 

In the above soft breaking terms we have omitted the family indices / and J. The above scalar fields 
correspond to the superfields which were indicated in our notation by a "hat". A^, A^ (with a = 1, 2 
or 3) and A^. (with 6 = 1, ... ,8) are the gauge superpartners of the U(l), SU(2) and SU(3) gauge 
bosons, respectively. 

The conventional wisdom is to assume a complete universality of the soft supersymmetric pa- 
rameters at the GUT scale. That is, a common scalar mass mo, a common gaugino mass M1/2 and 
a universal boundary condition for the trilinear soft breaking terms Au = Ad = Ae = AP . This 
gives rise to an appealing theory with a minimum number of free parameters, which, among other 
interesting features, allows Radiative EW Symmetry Breaking (REWSB) [Q. Such a SUGRA model 
has only two independent CP- violating phases at the GUT scale |T3| . These two physical phases can 
be chosen as arg(A'^) and arg(/i'^). Nonetheless, although very attractive, a complete universality 
of the soft breaking terms at the GUT scale (motivated by a "would be" GUT-scale SUGRA model 
scenario) is a simplifying assumption. If the universality of the trilinear soft breaking terms men- 
tioned above is relaxed, REWSB can still occur. Yet, in this less constrained framework, the At,Au 
and Ad (the only three soft trilinear breaking terms relevant for our discussion in this paper) can 
have different phases at the GUT scale, and therefore their phases at the EW scale are completely 
undetermined. As a consequence, we will show below that, contrary to previous claims, with no "fine 
tuning" of a.Tg{Au) and arglAd), the NEDM can meet its experimental limit even with squark masses 
of a few hundred GeV, leaving aig{At) large enough to drive significant CP-violating effects in top 
quark systems. We will discuss in detail the "SUSY CP problem" of the NEDM in section 4. 

Till now, no supersymmetric particle has been discovered, and in spite of the very fascinating 
theoretical features of the GUT-scale SUGRA model, there is no experimental hint that can indicate 
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the real nature of the underlying, short distance SUSY model. We therefore wish to keep an open 
mind and instead of restricting ourselves to the GUT-scale SUGRA model with complete universality 
of the soft breaking terms, our strategy will be to choose a plausible set of the SUSY parameters 



at the EW scale subject to the present experimental limits on the sparticle masses |16]. Indeed 



by relaxing the GUT-scale universality of only the soft breaking trilinear A terms, almost any low 
energy SUSY spectrum can be consistently recovered. 

We now describe the key features and assumptions of our low energy SUSY scenario: 
I) Motivated by the strong theoretical prediction of the unification of the SU(3), SU(2) and U(l) gauge 
couplings when SUSY particles with a mass scale around 1 TeV are folded into the Renormalization 
Group Equations (RGE), we will follow only one traditional simplifying GUT assumption; that there 
is an underlying grand unification. This leads us to have a common gaugino mass parameter defined 
at the GUT scale. Then, the difference between the three low energy gaugino mass parameters comes 
from the fact that they undergo a different renormalization when they evolve from the GUT scale to 
the EW scale due to the different gauge structure of their interactions. The gaugino masses, at the 
EW-scale, are then related by fl% : 



= sm 9w — = — , (7) 

a a as 

where 6w is the weak mixing angle and ma is the low energy gluino mass (from now on, we will refer 
to the gluino mass as mc)- Thus, once the gluino mass is set at the EW scale, the SU(2) and U(l) 
gaugino masses 1712 and mi, respectively, are determined. 

II) All the squarks except for the light stop are assumed to be degenerate with a mass Ms, while the 



light stop mass, which we denote by mi, is chosen subject to mi > 50 GeV ||T6[| . 

III) The mass matrices of the neutralinos, M^o, and charginos, M^, depend on the low energy Higgs 
mass parameter /i, the two gaugino masses m2 and rhi (which are resolved by the gluino mass) and 
tan/3 (see appendix A). Therefore, once /x, mo and tan/3 are set to their low energy values, the four 
physical neutralino species m^o [n = 1 — 4) and the two physical chargino species m^^ (m = 1, 2) 
are extracted by diagonalizing M^o and M^. In appendix A we list the full analytical prescription 
for diagonalizing the neutralinos and charginos mixing matrices in the case that fi is real. 

IV) We limit the lightest neutralino mass to be above 20 GeV and the lighter chargino mass to be 
above 65 GeV [g. 

V) We will always choose m^ > 200 GeV within the range M, - 200 GeV < m^ < + 200 GeV 
as roughly indicated by recent results from the Tevatron WE . 
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With no further assumptions, the low energy SUSY mass spectrum is fully determined from the 
only five free parameters /i, tug, M^, mi and tan (3. The possible CP phases, that can arise in this low 
energy SUSY framework, will be described in detail in section 4. 



3. Derivation of the CP Violating Asymmetry 

We now turn to the discussion of the possible CP-violating effects in the process pp ^ tb + X 
W^bb + X (i.e., production of a pair of tb followed by the top decay t — > bW^), driven by SUSY 
CP- violating sources |18[. To estimate the possible size of the effect, we consider the following 
CP-violating quantity: 

^ a(s) - a(s) 



a{s) + a{s) 

where s = {pt +^5)^ and a{s) and a{s) are the cross- sect ions for the reactions ud ^ tb ^ W^bb and 



ud ^ tb ^ W bb, respectively. Folding in the parton densities in the usual manner [|T^] we define 
the asymmetry of cross-sections at the level of the colliding protons as: 

a{sj + a[s) 

Now s = {pp + ppY and a{s) and a{s) are the cross sections for pp ^ tb + X ^ W^hb + X and 
pp ib + X W~bb + X, respectively. 

In the narrow width approximation the tree-level cross-section at the quark level is given by: 

/4sm Uws [s — my^^)^ 

In Fig. 1 we have drawn the tree level diagram for the parton level reaction ud — » tb. The possible 
super symmetric one loop diagrams that can violate CP invariance in this reaction are depicted in 
Fig. 2. Note that in both Fig. 1 and 2 the subsequent top decay is not explicitly drawn. The one-loop 
CP-violating triangle diagrams a-d of Fig. 2 in the production vertex will enter also the top decay 
t W^b with the t and momenta reversed. 

To one-loop order in perturbation theory, where the CP-violating virtual corrections enter only 
either the production or the decay vertices of the top in the overall reaction pp —>■ tb-\-X —>■ W^bb+X , 
and in the narrow width approximation for the decaying top, the overall CP asymmetry. A, can be 
broken into: 
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A = Ap + Ad. (11) 

In this Eq. Ap and Aj:, are the CP asymmetries emanating from the production and decay of the 
top, respectively, and are defined by: 



^ _ o{pp ^fb + X)- cr{pp ^ ib + X) 
^ ~ a{pp ^tb + X) + a{pp ^ib + X) ' ^ ^ 

_ Tit^Wn)-T{i^W-b) 
D Y{t^W+b)+T{t^W-b) ' ^ ' 

The PRA Ad^ defined in Eq. is an independent quantity and does not depend on the specific 
production mechanism of the top. We will therefore treat the CP-nonconserving asymmetries Ap and 
Ap) separately keeping in mind that the total CP-odd effect is simply the sum of the two. Moreover, 
it is convenient to further divide the cross-section asymmetry in tb production, Ap, into: 

Ap = + , (14) 

where X^^'^^^'^^ and A^p^ are the CP- violating cross-section asymmetries which arise from the triangle 
diagrams a — d and a' — d' and box diagrams e — h in Fig. 2, respectively. Similarly, the corresponding 
parton level asymmetries will be denoted by a "hat" . 

Asymmetry From Triangle Diagrams In The Production Amplitude - ^p'^'^sie. 

The one-loop tb and ib currents of the production amplitude can be parameterized as: 

y P=L,R \ J 
jli{tb) .gw _ f'Pl{k)Pb , a\ r> 

p¥p/' \^ + ^^^^^ W ' ' ^ ^ 

where P = L or R and L{R) = (1 — (+)75)/2. The index k indicates the triangle diagram (i.e. 



k = a,b, c, d, a', b', d or d'). V\^^ and V^i^-^ defined in Eqs. ^ and |T6|, contain the SUSY CP- violating 
phases as well as the absorptive phases of the triangle diagrams, both of which are needed to render 
^triangle ^Qj^.^gj-Q. Then, in terms of the scalar ('Pi(fc)) and vector {V2(k)) form factors, the parton 
level cross-section asymmetry ^^^^'^^^'^ jg given by: 
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i triangle 
Zip 



(17) 



where x = mf/s and the sum is carried out over all triangle diagrams in Fig. 2. It is easy to show 
that if one defines: 



J5. 



i(fc) 



X e 



iS 



i(fc) 



~ e 



iS'. 



2(fc) 



X e 



2(fe) 



(18) 
(19) 



where Sl^''\ d^^''^ are the CP-even absorptive phases (i.e., final state interaction phases) and Slj''\ 5'^^^ 
are the CP-odd phases associated with each of the triangle diagrams in Fig. 2, then: 



I i(k) 



X e 



-iS, 



i(fe) 



2(fe) 



X e 



-iS 



2(fe) 



(20) 
(21) 



The CP-violating scalar form factors in Eq. O, for each triangle diagram are then given by: 



Re(Pf(„)+P; 



1(a) 



Re(Pf(,)+P^J,)) 



1(c)) 



1(d), 



^^m,m,0^(Co" + Cr,) 

O TT 



a 



TT sin 6'vi/ 

a 

TT sin^ 



mi 



mAi(cr2 + c2^3) + ^x™a'c: 



1(c) J 



2ric 

12 



(22) 
(23) 

(24) 
Cj) ,(25) 



and the CP-violating vector form factors in Eq. O are: 



^<Ka) 








Re(P2'(.) 









a 



■ 2 n ^b^2A 1 
71 Sm U]y 



a 



w 



■5^22 



2 vr sin^ 9 



12 



(26) 
(27) 
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-^2K)) = 






Re(P2'(c') 


-^2^)) = 




-^2^)) = 



Re(Pi(,)-Pi(c))(^x?. 



Am ' ATI 



TTsin d-w 
1 a 



2 TT sin^ ^vi' 



Xm5 '^X" 



^X^^^c' 



(28) 

..y ' (29) 
(30) 

(31) 
(32) 



The SUSY CP-weak phases associated with the triangle diagrams, the O^'s above, are given in 
appendix B and all the corresponding loop form factors Cq , C^^ {p = 1,2 and g = 1 — 4) are defined 
and given in appendix C. 

Note that for diagrams a',b', d and d! only a vector form factor in the ud one-loop current 
contributes to CP- violation in the limit m^, (in this limit diagram a! contain no CP- violating 

phase). That is: 



_ .gw 



ij.(ud) _ .gw 



P=L,R 

y ^ p=j,.R 



(34) 
(35) 



and we can therefore consider these diagrams to have the same structure as diagrams a, 6, c and d 
(see Eqs. |l^ and |16D, taking the vector form factor in the ud current to be the vector form factor in 
the th current. 



Asymmetry From Box Diagrams In The Production Amplitude - Ap°^: 



There are two kinds of amplitudes for the box diagrams. The amplitudes for diagrams e and / 
in Fig. 2 can be written schematically as: 



whereas for Diagrams g and h: 



Vb 



(36) 
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-iM 



k=g,h 



167r2 



/13 /14 



■ (37) 



The four vectors pi and "pi-, that appear in Eqs. ^ and ^ above, stand for either q - the integrated 
four momentum in the loop or p = + pf, - the s-channel four momentum. Also, the denominators 
JAy. of the particles in the loop can be read off appendix C (see also pO|). 

A straight forward calculation of the interference of diagrams e, /, g and h with the tree level 
diagram yields: 



A 



box 



2 7rsin^ Bw {\-x)\2^x\ 



1 /•t+=o 

Ell 



t~=mi—s 



(3^ 



where y = m^^/s and the three Mandelstam variables at the parton level (i.e., s, i and u) are related 
via (for nib = = = 0)s + t + u = m^. For later use, we furthermore define the following 
quantity: 



X = - 
2 



{s - mf) - i (i - mf^ + n (m - )] . (39) 

For diagrams e and /, in the limit rriu = rria = m^, ^ 0, the X^, X^^, y^, and 3^^^ terms do not 
contribute to and we get: 



i Tr [LMiMi - 



whereas, for diagrams g and /i the , Xl/", y^, y'j^ and y'f^ terms vanish, and we have: 



(40) 
(41) 



-2mtulm {Xl^'^y^) Im ' ^ "^^"^ 



Zll Zl2 A3 M 
•^k-^k-^k-^k , 

Tr \Ljij)a 



1 fc'4.^fc/4.^fc/4.^t^4.^ 



(42) 
(43) 
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Blis,i) 



B}{s,i) 



-2Xmtm^<^Ol{Dl + Dl,) , 

20l [XmlDl^ + [u{u - ml){s - ml) - mlui 

+ X{u - ml)Dl^ + 2u{u - ml)Dl^ 

mluDl^ + u{s - ml)Dl2 + XD^^ 



^2A 



^2 [(X ■ 



B'Mi) ( 
0, 



u[u 



m. 



D 



11 



An 



MM 



ol 



(44) 



(45) 
(46) 
(47) 
(48) 
(49) 
(50) 
(51) 



The SUSY CP-weak phases associated with the box diagrams, the 0fc's above, are also given in 
appendix B and the corresponding four point one-loop form factors Dq, D^^ {p = 1,2 and g = 1 — 7) 
are defined and given in appendix C. 



PRA From The Decay Amplitude - A^: 



As mentioned before, there are four potential diagrams that can give rise to CP-violating PRA 
in the decay t W^h within the MSSM. These are diagrams a-d in Fig. 2 where it should be 
understood that the t and momenta are reversed and the ud current is disregarded. It is 
convenient to denote these one- loop top decay diagrams with the same alphabetical order, i.e., a, 6, c 
and (i, similar to the corresponding production diagrams, as the SUSY CP-odd phases defined in 
Eqs. p2| - p9| , the {, ^ ^^'s, are the same for both production and decay of the top. 

We parameterize the t W^h and i W~b decay vertices as follows: 



V ^ P=L,R 



1 1 



V P=L,R 



mt 



(52) 
(53) 
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Similar to the production process discussed above, I'^'^ and V^^^^ defined in Eqs. |5^ and contain 
the CP- violating odd phase as well as the absorptive phases of the decay diagram k, k = a — d. Then, 
in terms of the scalar {T>i(^k)) and vector (T'2(A:)) form factors, the PRA, Ajj, for the decay process is 
given by: 



A 



D 



E 



[X 
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Re(Pf(,)+PS,))+Re(D^ 



V. 



2(fc)J 



(54) 



k ^2(x + 2)- 

Now X = m1/m^ and the sum is carried out over all decay diagrams in Fig. 2 (i.e., k = a,b,c and 
d). The relations between and '^^k) between 1^2{k) ^2{fc) same as the relation 

between P^^) and V^f.^ and between 'P2(A:) '^2{k)y respectively, as is given in Eqs. 
For the scalar form factors in Eq. p3 we get: 



Re(Pf(„)+Pf(,)) 
Re(Pf(,)+:Df(,)) 
Re(Pf(,)+:Df(,)) 



Re(Df(,) + Vi 



1(d). 



3 vr 



a 



vrsin 9w 

a 



ITLt 



12 



-mt 



71 sin^ 9w 
Re(Df(,. + 



mtOl {C: 



1(c); V"^x° 



(55) 
(56) 

(57) 
Ci) ,(58) 



while the vector form factors in Eq. are given by: 





- '^2{a) 




- ^2(fe) 




- ^2(0 


Re(P2^(,) 


- ^2(d) 







24 5 



("^t (C'22 " C'23) + "^w^ (C'li + C'22 — C'12 ~ C'aa) + 2C. 



vr sin^ Oy/ 

1 a 

2 vr sin^ Oy/ 

+mtm^^Ol {CI, - CI,) - m,m^Ol (Q + CI, - CI, 
m^oOtCl 

2(c)) (' 



— m 



An 



ATI 



(59) 
(60) 



24 



C 



(61) 
(62) 



As mentioned above, the SUSY CP- weak phases for the decay diagrams, the O^'s above, are the same 
as those for the production triangle diagrams a-d in Fig. 2 and are given in appendix B. Although 
the same notation (as in the production case) for the corresponding loop form factors Cq, C^^ is used 
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in the decay case at hand, these form factors are separately defined and given in appendix C. 
4. CP Phases of The Low Energy MSSM and the Neutron EDM 



Before presenting our numerical results, we wish to clarify our approach with regard to the CP- 
violating sector of the low energy MSSM and their effect on the NEDM. Disregarding for now the 
possible phases at the GUT scale, we assume for simplicity that at the EW scale, all CP-violation 
(apart from the usual SM phases) resides in the complex trilinear soft breaking terms, the A/'s. 
For definiteness, we assume arg(/i) = motivated by analyzing the NEDM, which suggests that 
arg(/i) < O{10~^) is required for the NEDM to satisfy the experimental bound for squark masses of 
a few hundred GeV (which we are assuming throughout the paper) 0, ||, 0. 

In fact, from a phenomeno logical point of view, this is certainly a plausible scenario that can 
emanate from a GUT-scale SUSY model in which the universality of the A terms is relaxed to give 
arbitrary phases aTg{A^) at the GUT-scale. In this regard, we take note of a recent very interesting 
work by Garisto and Wells |^. They obtained severe constraints on the low energy phases of At 
and fi by deriving relations between the two. For that they use the complete set of RGE involving 
the complex parameters of a GUT-scale SUSY model with and without universal A terms and 
with some definite boundary conditions at the GUT-scale (for example arg(/i'^) = aig{A^) = 0). 
However, taking non-zero phases at the GUT-scale, the constraints obtained in may not hold 
as ,in general, assuming arbitrary phases at the GUT-scale one can practically make no prediction 
on the corresponding phases at the EW-scale. Our approach, the EW —>■ GUT approach, will be 
to assume a set of SUSY phases at the EW-scale, subject to existing experimental data, which 
implicitly assumes arbitrary phases at the scale in which the soft breaking terms are generated. 
With only the low energy phases of the various Af terms, if all the squark masses except for the 
stops are degenerate with a mass Ms, then only arg(At) contributes to CP- violation in the reaction 
pp ^ tb + X ^ W^bb + X. Any CP-odd phase from the trilinear soft breaking terms of the other 
squarks will enter into the asymmetry A, defined in Eq. ^, with a suppression factor of m^/m^ [q 
stands for any quark but the top) even if the squarks masses were not taken exactly degenerate. 
One power of niq/mt comes from the Dirac algebra when evaluating the squared matrix elements, 
if ntq is not neglected. Another irtq/mt comes from the squarks mixing matrices: as can be seen 
from appendix A. The CP-violating quantity that arises from qi — q^ mixing (i.e., Im(Z^**Zg*)) is 
proportional to the SM quark masses, thus for rriq/Aq —>■ there will be no mixing in the sector of 
the super symmetric partners of the light quarks. 
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In appendix A we list the Feynman rules needed for calculating the CP-odd effect in the reaction 
ud ^ tb ^ W~^bb. It is then obvious, that in principle, all the CP-violating vertices of the MSSM 
arise from diagonalization of the complex mass matrices Z^, ^at, Z~ and Z~^. However, with 
the assumption of arg(/i) = 0, it turns out that Zn, Z~ and Z'^ are real. In particular, assuming 
the universality of the gaugino masses at the GUT-scale, the phase of the common gaugino mass, 
arg(Mi/2), can be set to zero by a phase rotation |^, thus leaving the fhi i = 1 — 3 "phaseless" 



at any scale. We therefore have all CP-violation arising from — mixing (i.e., from aig{At))\ in 
particular from (see appendix A): 

Im(ei) = \m{ZY*Zf) = '^"^ '^"^ • 

We will choose maximal CP-violation in the sense that \m.{Q) = (— 1)*~^/2 thus presenting all our 
numerical results in units of sin2^t sin jSt- With no further assumptions, what is left to be considered 
is the NEDM for a chosen set of the free parameters of the above phenomenological low energy 
MSSM. 

The EDM of the neutron, d^, is presumably one of the most important phenomenological problems 
associated with SUSY models especially with regard to CP violation. With a low energy MSSM that 
originates from a GUT-scale SUGRA model (with complete universality of the soft breaking terms), 
keeping dn within its allowed experimental value (i.e., dn < 1.1 x 10~^^ e-cm [Q) requires the "fine 
tuning" of the SUSY phases to be less then or of the order of 10^^ — 10^^ for SUSY particle masses 
at around several hundred GeV's 0. 

When arg(//) = 0, the leading contribution to a light quark EDM comes from gluino exchange, 
which, with the approximation of degenerate u and d squark masses (which we will denote by rrig), 
can be written as 0, 0] : 

, 2q!s ^ |y4o|sinao . , 

where mglmq) is the quark(squark) mass and Qg is its charge. Also, r = rriQ/m'^ and K{r) is given 
by 

K(r) = ^^-^(- + -r +^-^^±^\nr] . (65) 
^ ^ (r- 1)3 \^2 2 1 - r y ^ ' 

Aq is the complex trilinear soft breaking term at the EW scale associated with the squark g, and we 

have defined Aq = lAgle*"'. Then, within the naive Quark Model, the NEDM can be obtained by 

relating it to the u and d quarks EDMs (i.e., du and dd, respectively) by dn = {Md — du)/3. 
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We stress again that with no universal boundary conditions for the soft trihnear terms and their 
phases at the GUT scale (i.e., arg(Af/) = aig^Ai)) = arg(A£;) = arg(A'^)), there is no a-priori 
reason to believe that the low energy phases associated with the different Af soft breaking terms 
are related at the EW scale. We therefore consider arg(y4„) and arg(Ad) to be free parameters of 
the model no matter what aiglAt) is. In Figs. 3a and 3b we have plotted the allowed regions in 
the sina^j — sina^ plane for not to exceed 1.1 x 10~^^ e-cm (the present experimental limit) 
and 3 x 10~^^ e-cm, respectively. In calculating dn we assumed that the above naive Quark Model 
relation holds. Although there is no doubt that it can serve as a good approximation for an order of 
magnitude estimate it may still deviate from the true theoretical value which involves uncertainties 
in the calculation of the corresponding hadronic matrix elements (see and references therein). 
Note also that it was recently argued that the naive quark model overestimates the NEDM, as the 
strange quark may carry an appreciable fraction of the neutron spin which can partly screen the 
contributions to the NEDM coming from the u and the d quarks [0]. To be on the safe side, we 
therefore slightly relax the theoretical limit on c/„ in Fig. 3b to be 3 x 10~^^ e-cm. 

We have used, for these plots, fn^ = = Mg = 400 GeV, rriG = 500 GeV and for simplicity we 
also took \Au\ = \Aa\ = Mg [Mg is the only high SUSY mass scale associated with the squarks sector 
in our low energy MSSM. Therefore it is only natural to choose the mass scale of the soft breaking 
terms according to Mg). Also, we took the current quarks masses as = 10 MeV, m„ = 5 MeV 
and ag{mz) = 0.118. 

From Fig. 3a and in particular Fig. 3b, it is evident that Mg = 400 GeV and = 500 GeV 
can be safely assumed, leaving "enough room" in the sina„ — sina^ plane for \dn\ not to exceed 
1.1 — 3 X 10~^^ e-cm. In particular, while sina^j is basically not constrained, depending on sina^, 
—0.35 < sinorf < 0.35 is needed for \dn\ < 1.1 x 10~^^ e-cm and —0.55 < sina^^ < 0.55 is needed for 
\dn\ < 3 X 10~^^ e-cm. Moreover, varying mc between 250 GeV to 650 GeV (we will vary mo in this 
range when discussing the numerical results below) almost has no effect on the allowed areas in the 
sin a„ — sin plane that are shown in Figs. 3a and 3b. That is, keeping Mg = 400 GeV and lowering 
mo down to 250 GeV, very slightly shrinks the dark areas in Figs. 3a and 3b, whereas, increasing 
niG up to 650 GeV slightly widens them. Of course, dn strongly depends on the scalar mass Mg - 
increasing Mg enlarges the allowed regions in Figs. 3a and 3b as expected from Eq. |6^. It is also 
very interesting to note from Fig. 3a that, in some instances, for a cancellation between the u and 
d quarks contributions to apply, sina„,sinarf > 0.1 is essential rather than being just possible. For 
example, with | sinaul > 0.75, | sina^l > 0.1 is required in order to keep dn below its experimental 
limit. 
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We can therefore conclude that CP-odd phases in the A terms at the order of few x 10~^ can be 
accommodated without too much difficulty with the existing experimental constraint on the NEDM 
even for typical SUSY masses of < 500 GeV. Therefore, somewhat in contrast to the commonly 
held viewpoint we do not find that a "fine-tuning" at the level of 10~^ is necessarily required for the 
SUSY CP-phases nor for the squark masses. 



5. Numerical Results 

We now turn to the discussion of our main numerical results. We first focus on CP-violating 
effects in the production amplitude which gives rise to the cross-section asymmetry, Ap, and then we 
will present our results for the PRA, Ad, associated with the top decay t — > W~^b and its conjugate 
one. Instead of using the asymmetry Ap defined in Eq. after folding in the parton luminosities in 
the usual manner (see [|l^]), we define a Partially Integrated Cross-section Asymmetry (PICA) with 



respect to the variable r or equivalently to s (s = ts). Thus: 



^PICA _ 



dr 



a{s = Ts) — a{s = Ts] 



^^^^ [a{s = TS} + a{s = TS 



dT 

— , (66) 

CLT 



(^a{s = ts) + a{s 

where we have introduced the "parton luminosity" dCud{T) / dT ||19|. We also define the invariant mass 
of the tb system to be rritb = Then, in the numerical evaluation of the asymmetry A^*^^ we choose 
r_ and r+ according to = rrit {nib = is assumed through out the paper) and = 350 GeV, 
respectively. In the Tevatron, one of the problematic "backgrounds" to the reaction pp ^ tb + X is ti 



production |23]. In principle, the tt background can be eliminated by imposing the above naive upper 
cut on rritb, which is below the ti production threshold (i.e., = 350 GeV < 2mt). Although the 
actual experimental cuts that will be made in order to remove the tt "background" in the Tevatron 
(when studying the tb final state) might be more involved, this naive cut on rritf, serves our purpose. 
A detailed discussion of the exact experimental cuts for the tb final state is beyond the scope of this 
paper but can be found in [|I], Note also that this upper cut on rritb has practically no effect on 
the cross-section for pp —>■ tb + X (or equivalently the available number of tb events in the future 
2 TeV Tevatron), as most of the tb pairs will be produced with an invariant mass of rritb S 350 
GeV. In particular, the fully integrated cross-section for pp tb + X is (i.e., up to rritb = 2 TeV) 
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a{s = 2 TeV) ~ 360 [fb] out of which ~ 80% of the tb pairs are produced with an invariant mass 
mt, < 350 GeV. 

For the numerical evaluation of the CP- violating asymmetry A^'^^, we take Ms = 400 GeV. 
We also vary the gluino mass between 250 GeV < tug < 650 GeV (as was emphasized in the 
previous section, these masses of the degenerate squarks and the gluino do not contradict the existing 
experimental upper limit on the NEDM). We choose two representative values for tan (3; tan (3 = 1.5 
and tan (3 = 35 which correspond to a low and a high tan (3 scenario, respectively. The low energy 
Higgs mass parameter fi is varied between -400 to +400 GeV, whereas, the mass of the light stop 
(mi) is varied between 50 to 400 GeV. As mentioned above, we also choose maximal CP-violation, 
driven by t/, — tR mixing, in the sense that Q = (— 1)^~Y2 (see Eq. |63|), whereas, we take arg(yu) = 0. 
Therefore the asymmetries are always given in units of sin 26t sin jSt- 

The consequences of this low energy MSSM scenario (with the above chosen mass spectrum and 
CP-odd phases) with regard to the various diagrams that are depicted in Fig. 2 are: 

I) The triangle diagrams 2c, 2b', 2c' and 2d' as well as the box diagrams 2/ and 2h do not acquire 
any CP-violating phase in the limit arg(/i) = 0. 

II) Diagram 2a' does not have any CP-violating phase in the limit 771^,171^ —>■ 0. 

III) Diagram 2a does not have an absorptive cut for rritb < 350 GeV. 

IV) Diagram 2b is p-wave suppressed near threshold. Thus its contribution to A^^^ is negligible 
compared to the other diagrams. 

With the above points I-IV, the study of CP- violation in this reaction simplifies to a large extent. 
That is, we are left with only three CP-violating diagrams that can make important contributions 
to A^'-^^. These are: diagram 2d, 2e and 2g out of which diagram 2d is expected to be the most 
important one. In particular, in the range where the asymmetry exceeds the 1% level, diagram 2d 
is responsible for ~ 90% of the total CP-odd effect in Ap^*-^"^. Moreover, it is important to note that 
while the CP-violating effects from diagrams 2e and 2g strongly depend on the heavy scalars mass 
Ms (through two scalar exchanges in the loop), diagram 2d is less sensitive to Mg having only one 
scalar exchange in the loop. In particular, as Mg increases, diagram 2d becomes relatively more 
dominant as far as CP is concerned (Of course, in general, increasing the SUSY scale Ms causes the 
asymmetry to drop for diagrams 2d as well as for diagrams 2e and 2g). It is worth mentioning that, 
in that range where the asymmetry is above the percent level mainly due to the CP-effects coming 
from diagram 2d, the CP-violating effect driven by diagrams 2e and 2g has the same relative sign to 
that of diagram 2d, thus yielding a slightly bigger overall asymmetry. Nonetheless, for simplicity, we 
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will present all our numerical results only for diagram 2d. In any case, we are only interested in a 
rough estimate of the CP- violating effect in the reaction pp ^ tb + X and not in the exact numbers, 
as many simplifying assumptions had to be made along the way. 

In Figs. 4a and 4b we have plotted Ap^"^ for Mg = 400 GeV, mi = 50 GeV and for three 
values of itlg as a function of the low energy Higgs mass parameter fi, for tan (3 =1.5 and tan/3 = 
35, respectively |2^, We see from Fig. 4a that in the low tan/5 scenario and a gluino mass 
between 350 GeV < uig < 550 GeV, moderately negative values of fi, lying in a narrow range, 
—100 GeV < /i < — 70 GeV, are required for the asymmetry to be above 2%. In the best case, 
for tug = 450 GeV the asymmetry peaks at around fi ~ —90 GeV, reaching ~ 2.75%. A positive /i 
around 140(240) GeV may also give rise to a ~ 2(1)% asymmetry for a gluino mass niG ~ 550 GeV. 
Also, for fJ' < — 200 GeV and fi > 350 GeV the asymmetry drops below the 0.25% level. From 
Fig. 4b (high tan(/?) we see that around ~ 110 GeV the asymmetry can also reach the 2% level 
for 400 GeV < mc < 575 GeV (see also Fig. 6b). Here, the asymmetry drops below the 0.25% level 
for l/^l ^ 250 GeV. Notice that, in the high tan/5 scenario, A^'~^^ is almost insensitive to the sign of 
/i. This happens mainly due to the fact that, for high values of tan/5, the charginos and neutralinos 
masses are almost independent of the sign of fi as can be also seen in Table 1. The only terms in 
the chargino and neutralino masses which linearly depend on /i are proportional to sin 2/5 which is 
of the order of few x 10"^ for tan/5 > 30 (see Eq. A.22 and A.25-A.32). 

Before discussing the results in Figs. 5-11 we remark that, in what follows, we choose several 
representative values of /i; some of them maximize CP-effect in the production amplitude and some 
maximize the CP-violating PRA effect in the decay t — > bW to be discussed later. In particular, for 
tan/5 = 1.5 we always take the values fi = —70, —90, —130, 140 and 240 GeV, while for tan/5 = 35 
we choose /x = -110,-170,-190,110 and 170 GeV, although, in some cases the combination of 
{tan/5, /i, mc} for some of those values is forbidden according to our criteria (see section 2). It is 
useful to keep track of the charginos and neutralinos masses for a given set of {tan /?, /i, mc}, especially 
due to the various absorptive thresholds that can emanate from the loop integral of diagram 2d |2^ . 
Therefore, for the reader's convenience, we give in Table 1 the masses of the charginos and neutralinos 
for various sets of {tan /?, /i, uig} that are repeatedly being used throughout this analysis. 

In Figs. 5a and 5b we show the dependence of A^'-^^ on mi, the mass of the light stop, for various 
values of ^, Mg = 400 GeV and for tan/5 = 1.5 and tan/5 = 35, respectively. However, the general 
behavior that is depicted in Figs. 5a and 5b holds for any value of fi. That is, as expected, when mi 
is increased the asymmetry falls till it totally vanishes for mi = 400 GeV in which case the two stop 
species are degenerate. It is evident from Figs. 5a and 5b that, in general, mi < 75 GeV is needed 
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for the asymmetry to be above ~ 1%, although in some cases, for example when // = —70 GeV and 
tan/? = 1.5, a 1% asymmetry can arise even with mi > 100 GeV. 

The dependence of A^^^ on the gluino mass rriG is shown in Figs. 6a and 6b for tan/? = 1.5 
and tan /? = 35, respectively. Here again, we took Ms = 400 GeV and nii = 50 GeV and choose the 
same representative values for /i. We see that for a small tan/? and a negative fi around —100 GeV, 
^P^^ > 2.5% for 350 GeV < < 550 GeV and peaks at around nic ~ 400 GeV. Note that even 
a heavy gluino, e.g., ~ 650 GeV, can give rise to a > 2% asymmetry if fi ^ —90 GeV and tan/? 
is of order one. We see from Fig. 6b that for large tan/? and for \fi\ = 110 GeV, A^'^^ ~ 2% for 
400 GeV < rriG ^ 550 GeV and peaks at itlg ~ 450 GeV. In both the low and high tan /? scenarios, 
with a "light" gluino, itlg < 300 GeV, the cross-section asymmetry is below 1%. 

Figs. 7a and 7b show the dependence of A^^^ on tan /? for rriG = 350 and 550 GeV, respectively. 
Evidently, as far as CP is concerned, a low tan /? is better then a high one. We see that as one goes 
to tan/? > 10, the asymmetry is almost insensitive to tan/?. This also holds for high values up to 
tan /? = 65 which are not shown in Figs. 7a and 7b. 

Until now we were not interested in the overall sign of the cross-section asymmetry A^^^. Includ- 
ing the PRA CP-violating effect, A^, coming from the decaying top, the relative sign between A^^^ 
and Ao becomes important as the total CP-violating effect is the sum of the two i.e., A = A^^'^^ + Ad- 
However, we will show below that, in most instances, Ao is smaller than Ap^*"^ by more than an 
order of magnitude and, therefore, its relevance to the CP-effect in pp tb + X ^ W^bb + X is 
negligible. 

Let us now discuss the PRA effect, A^. Once again, we will use our low energy MSSM scenario 
(described in detail in section 2 and 4) to estimate the CP-odd effect. Although, as was mentioned 
before, the PRA Ad does not depend on the specific production mechanisms of the t and t, in most 
cases we will evaluate its magnitude within the same ranges of the SUSY free parameter space as was 
taken in Figs. 4-7. This should enable the reader to easily extract the overall asymmetry emanating 
from both production and decay processes discussed in this paper. However, as was mentioned 
above, throughout most of the range of our SUSY parameter space, A^ is more than one order 
of magnitude smaller than A^p*^^, therefore, the relative sign between A^ and A^p*^^ is essentially 
irrelevant. Note that the consequences of the low energy MSSM framework on the various diagrams 
(i.e., the "reversed" diagrams a — d in Fig. 2) that can potentially contribute to Ao are: 

I) For niG > 250 GeV diagram 2a does not have the needed absorptive cut and thus, does not 
contribute to Ao- 

II) As in the case of production, diagram 2c does not have a CP-violating phase in the limit arg(/x) = 
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0. We are therefore left with only two diagrams that can contribute to A^- These are: diagram 2b 
and 2d, where the leading contribution to again comes from diagram 2d. 

Our main results for A^ are shown in Figs. 8-11. A quick look at Figs. 8-11 reveals that, 
throughout a large portion of the SUSY parameter space discussed here, Ad tends to be below 0.1%, 
although for some specific values of the parameters it can approach 0.3%. The asymmetry we find 
is therefore somewhat small compared to the estimates of Grzadkowski and Keung (GK) and of 
Christova and Fabbrichesi (CF) in 0. In the GK limit only the gluino exchange of diagram 2a 
was considered in which case < 120 GeV is required in order to have the necessary absorptive 
cut (when mi = 50 GeV). In the best case, GK found a ~ 1% asymmetry for = rrii = 100 
GeV |26[. On the other hand, in the CF limit, numerical results were given only for the neutralino 
exchange diagram (i.e., diagram 2b) wherein the CP-phase was chosen to be proportional to arg(yu) 
and maximal CP- violation with regard to arg(/i) was taken [^. In the best case, CF found a ~ 2% 
asymmetry for = 100 GeV. However, each of those largish PRA asymmetries, reported by GK 
and CF in 0, suffer from two drawbacks, (i) For the GK limit, mc < 120 GeV is now essentially 
disallowed by the current experimental bounds, (ii) For the CF limit, arg(/i) > 10~^ is an unnatural 
choice bearing the stringent constraints on this phase coming from the experimental limits on the 
NEDM as discussed in section 4. {iii) For both the GK and CF limits, the large asymmetry arises 
once the masses of the superpartners of the light quarks are set to 100 GeV. Again, this is a rather 
unnatural choice as it is theoretically very hard, if at all possible, to meet the NEDM experimental 
limits when the masses of the squarks (except for the lighter stop) are at the order of 100 GeV. 
Instead, we discuss below the PRA effect, A^, including all possible diagrams and subject to the 
available experimental bounds on both the NEDM and the supersymmetric spectrum. In particular, 
we again take arg(/i) = 0, rrig = = 400 GeV, rrii > 50 GeV, > 250 GeV, the LSP mass to be 
above 20 GeV and the mass of the lighter chargino to be above 65 GeV (see sections 2 and 4) [ p^ . 

As can be seen from Figs. 8a, 9a and 10a, in the low tan f3 case (i.e., tan (3 = 1.5) and for negative 
values of /i in the range —90 GeV < /i < — 30 GeV, Ad can be about 0.1-0.2% provided the gluino 
mass lies in the range 300 GeV < niG < 550 GeV and the lighter stop mass (mg) is between ~ 50 
and ~ 70 GeV. In fact. Ad can even reach ~ 0.3% if in addition to mi being in that narrow range we 
also have mc ~ 320 GeV. As fi becomes positive, at /i ~ 150 GeV, Ad tends to grow with the gluino 
mass. In particular, we find for example that for fi ^ 140 GeV and mo > 600 GeV, Ad ~ 0.2% 
becomes possible. The results for the high tan/5 case (i.e., tan/5 = 35) are shown in Figs. 8b, 9b and 
10b. Evidently, with a high tan/5, \Ad\ < 0.1% throughout all the range of our SUSY parameter 
space. 
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The dependence of on the hghter stop mass, mi, is shown in Figs. 9a and 9b. As expected, 
Ad drops as m; is increased and vanishes for mi > rrit — mLSP ('tilsp is the LSP mass), for which 
case there is no absorptive cut in any of the diagrams a — ci in Fig. 2 . 

Finally, in Figs. 11a and lib we evaluate Ad as a function of tan/5 for mc = 350 and 550 GeV, 
respectively. We can see that, here also, a small tan (3 gives rise to a bigger asymmetry and, as tan (3 
grows above ~ 10, Ad becomes insensitive to it. In the best cases, for tan/5 ~ \Ad\ ~ 0.3% 

become possible for both = 350 and 550 GeV. 

Before summarizing, we wish to emphasize that the decay asymmetry Ad is independent of the 
production mechanism. Therefore, the LHC with its large (~ 10^) production rate for ti pairs could 
be a good place to search for Ad- To probe such a PRA signal in the top decay t —>■ bW, of the order 
of a few tenths of a percent, will naively require ~ 10^ tt pairs when no efficiency factors are taken 
into account. Note that such a measurement requires only the detection of the charge of the top for 
which the systematic errors may be kept relatively small. Even if we take an efficiency overall factor 
of ~ 0.1, then with ~ 10^ ti pairs the reach of the LHC will be around Ad > 0.3%. Although, as we 
discussed in the preceding pages, in our SUSY study it appears difficult to attain that large a decay 
asymmetry, the experimental search would still be worthwhile. 
6. Summary and Conclusions 

To summarize, we found that within a portion of a plausible MSSM low energy parameter space, 
a CP-violating cross-section asymmetry in the reaction pp —>■ tb + X can be at the level of a few 
percent. Furthermore, a CP-violating PRA effect in the subsequent top decay cannot exceed the 0.3% 
level throughout our chosen range of the MSSM free parameter space. Therefore, the CP-violating 
asymmetry in this reaction arises predominantly from the production vertex. These asymmetries, 
in production and decay of the top, are driven by the complex entry of the soft trilinear breaking 
term associated with the top, aTg{At). In particular, we have shown that for scalar SUSY masses as 
well as a gluino mass of the order of 0.5 TeV, and with a Higgs mass parameter, /i, at the range of 
—250 GeV < /i < 250 GeV a CP-violating signal above the percent level might indeed arise in the 
reaction pp ^ tb + X —>■ W^bb + X, provided that the light stop particle have a mass below 100 
GeV. 

We have also shown that this phenomenologically acceptable low energy MSSM scenario does 
not contradict the experimental limit on the NEDM (i.e., (i„ < 1.1 x 10~^^ e-cm). In particular, in 
our low energy MSSM framework, SUSY CP-odd phases in the trilinear soft breaking terms (Af's) 
at the order of few x 10~^ are allowed even with squark masses below 500 GeV. The theoretical 
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uncertainties involved in calculating the NEDM may even strengthen this statement. In this respect 
our findings are somewhat different from what is commonly claimed in the literature. 

In the future 2 TeV pp Tevatron collider, the cross-section for pp tb + X is expected to be 
about 300 (fb) if a cut of rritb < 350 GeV is applied on the invariant mass of the tb. Therefore with 
an integrated luminosity of £ = 30 (fb)~^ |l|, |2^ an asymmetry ApicA ~ 3%, which was found here 
in the best cases, can be naively detected with a statistical significance of 3a. We must caution, 
however, that the values of the asymmetry as large as 3% appear attainable only for some specific 
range for some of the input parameters. Note also such a detection at the Tevatron will require the 
identification of all tb pairs , which, in principle, can be achieved only if the top can be reconstructed 
even when the W decays hadronically. Therefore, at the future 2 TeV pp collider, a percent level 
CP- violating signal in the reaction pp tb + X —>■ W~^bb + X may become accessible. We have also 
emphasized that the PRA effect in the decay t bW may be measurable at the LHC 

Note that the above line of reasoning can be reversed. That is, one can ask: what limit can be 
obtained on a.Tg{At) by studying CP-violation in the reaction pp ^ tb + X —>■ W~^bb + X at the next 
runs of the Tevatron or in the decay t bW at the LHC? We found that in the best case, the 2 
TeV Tevatron with C = 30 (fb)~^ will be marginally sensitive to CP-odd effects coming from arg(y4() 
and, taking into account the various efficiency factors for such a detection, perhaps a ~ 1 — 2-sigma 
bound on a.Tg{At) will be feasible. However, note that in a possible Tevatron upgrade with center 
of mass energy of 4 TeV the production rate of tb becomes three times bigger and, therefore, 
in such an upgraded version of the Tevatron a more stringent upper limit may indeed be placed on 
aig{At) through the study of the tb cross-section as suggested here. The LHC may also be sensitive 
to a.Tg{At) through a study of a PRA effect in the top decay t — > bW if the SUSY parameter space 
lies in the window wherein the PRA effect reaches its maximum values i.e., 0.1-0.3%. In such a 
scenario a statistically significant limit on arg(y4() may be within its reach. 

In closing, we wish to remark that it will be useful to explore the SUSY mediated CP-violating 
effects that can emanate in the l^-gluon fusion subprocess which contributes to the same final state 
(i.e., Wg — i> tbd) and which has a comparable production rate to that of the simple ud — > tb in the 
2 TeV Tevatron. Unlike in the case of the 2HDM where, to one-loop order, no CP- violating 
corrections enter the M^-gluon subprocess, within the MSSM, various one-loop triangle and box 
corrections can give rise to CP-nonconservation also in the H^-gluon fusion subprocess. However, 
it is also important to emphasize that due to the extra light jet in Wg — > tbd and the different 
kinematics, in principle, the tb pairs originating from the ly-gluon fusion can be distinguished from 
those produced via the ud fusion |2^ . 
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Appendix A 



In this appendix we list the relevant pieces of the SUSY lagrangian written in terms of the mass 
eigenstates. We also give the analytical expressions for all the corresponding mass mixing matrices 
Z/, Zi\i, Z+, and the masses of the supersymmetric particles which are functions of these matrix 
elements. 



The SUSY lagrangian pieces are |T^: 



djUXr, 



C- - 
r - 

'-'djdg 

^WXmXl 



^J_Z]^ Z^:V-\^ ut)W- + H.c. 



ufd 



V2 



u ^2m 



R + 



"^"^'zl^Z^^Llv-'rm + H.c, 



Vl 



J Am 



V2m. 

Vl 



d y2j rj- 

^d ^2m 



L- 



Z]^ZZR \ V'^'u + H.c. 



V2 
V'i Xn 



ud„ 



■'d 2m 



9 __2}i*L^ 



V2 " 



V2 



2i* ryAn 

N 



ZTZ 



L+ 



gi^ndwZ^'^ZT- 



2i* ryln* V^IT^U r7li* rvAn* 



An 



V2 



R}u + H.c. 



y2j y3n 

^d 



L+ 



-y^gtanewZ'.'Z'^-*- 



2j yln* _ V^ITT'd ylj y3n* 



Vl 



-Z^'*L + Z^"R 
-Z^L + z^M 



d 

U + H.C. 

d + H.c. , 



R\d + H.c. 



{K-L + K+R}x'iW; + H.c. 



{A.l] 



(A.2) 



(A.3) 



(A.4) 



(A.5) 

(A.6) 
(A.7) 
(A.8) 



where L{R) = |(1 — (+)75) and the Ui and u {dj and d) stand for up squark and up quark (down 
squark and down quark), respectively. Also Xm, Xn 9 charginos, neutralinos and gluinos 

respectively. We have also defined: 
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K- 



1 
3 



v3n* / 



1 



7Z(t 7+* 7+* 



and the mixing matrices Z^-, Z^, Z^^, Z and Z+ are given below. 

Let us define the following diagonalizing mass matrices ||^, |T2|, |13|, ^ 



(A.9) 
(A.IO) 
(A.ll) 



ZfMjZj - 
(Z-fM^Z^ 

ZjjM^oZjq = 



diag (mj-^ , m| 
= diag {m^,,m^^] 
diag ( 



m^o,m^o,m^o,m^o 



(A.12) 
(A.13) 
(A. 14) 



where Mj is the mass squared matrix of the scalar partners of a fermion. and M^o are the mass 
matrices of the charginos and neutralinos, respectively. 
Mj is then given by [0 0, |ll]: 



m) - cos 2(3 {Tsf - sin^ Ow)M^ + 
-mf{Rffi* + Aj) 



JL 



-mf{Rff^ + A}) 
ml - cos 2(3Q / sin^ OwM^ + m 



(A.15) 



where mj is the mass of the fermion f,Qf its electric charge and T^f the third component of the 



weak isospin of a left-handed fermion /. m| (m^ ) is the low energy mass squared parameter for the 



left(right) sfermion /L(/ij). Rf = cot/?(tan/?) for Ty = |( — |) where tan/? = V2/V1. 
and M^o are given by 0, 0, [T3|, [T^, 



m2 v^Miy sin /? 

V^MvK cos /3 /i 



mi 
m2 
-M^ cos /5 sin 6*1^ cos /5 cos 6'vi/ 

sin P sin 6*1^/ — sin /5 cos 6'iy 



— cos P sin 6*^1/ 
cos P cos 6*1^ 




Mz sin /5 sin 
-Mz sin /5 cos 

-/^ 




(A.16) 



(A.17) 



where mi (7712) is the mass parameter for the U{1){SU{2)) gaugino. 

After diagonalizing these mass matrices according to Eq. A.12-A.14, the corresponding masses 
of the supersymmetric partners of the fermions, the charginos and the neutralinos are given by: 
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+ 



b'j: + C 



/l{/2) 2 

df = 2m J — cos2/5T3jM| + + m^^ 
hf = cos 2/3(2Q/ sin^ Ow - T3f )Ml + m) 



2 

h 



m 



-2m t 



Rffi + A} 



(A.18) 
(A.19) 
(A.20) 
(A.21) 



+ 



[ml 



l?f + 4M^ cos' 2/5+ 



+ 4M^ (m^ 



/i^ + 2m2/usin2/?) 



(A.22) 



In Eq. A. 14 we have defined Zat such that the elements of the diagonal neutralinos mass matrix are 
real and non-negative. It is sometimes more convenient to allow for negative entries for the m^o 



which implies the change of Eq. A. 14 to 29, 30 : 



NM^^N = diag(eim^o, e2mj^o, 63^-^0, ^4771^0) , (A. 23) 

where is a real matrix. While the m^o are always positive, the e^'s are either ±1. With this 
substitution, in the Feynman rules one has to use the relation [^, ^ (note the slight difference 
between our notation for Zj^ and the one used in PD| ): 

^^.■ = (v^)*%. (A-24) 
The positive real masses of the neutralinos are then given by: 



ei"^x? = -«x + V^x + Cx + c^x ' (A-25) 



e2m^o = a^-^6--c^ + rf- , (A.26) 



egm^o = -a^ - ^h^ + + rf^ , (A.27) 



e4m^o = flx + V ^x ~ '^x + '^x ' (A. 28) 

where: 
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2 3 

C3 



= -(mi + m2) 



and a, C2 and C3 are defined in |]T^ with the substitution /x — ^ 
The mass mixing matrices are then given by 

cosOf —e~^^f sinOf 
e'^^^sinOf cos Of 

where: 



-/X. 



tan6'f = . 



When /i is chosen as real, the phase /?/ is given by: 

tan/3/ 

where = arg(yl/). 







sin 9aj: 


fiRf + \Af 


cos 9Af 



The charginos mixing matrices are given by 0: 

(0+)-i i/ det(M^) > 
(a30+)-i 2/ det(M^)<0 



o^, z- 



where: 



Oh 



cos 0-1- sm (p± 
— sin 0± cos (f)± 



and: 



tan 20_ = 2v/2Mi 



/i sin /? + m2 cos P 
ml-fi'^~ 2M^ cos 2p 
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H cos (3 + 1712 sin P 



The neutralinos mixing matrix elements are given by |[T7| : 



(A.39) 



i2 



Nil 



1 mi — eim^o 

Nil tan 9w ^2 — ^im^o 

Nis 
Nil 



(A.40) 



Ii[m2 — eim^o][mi — ejm^o] — M| sin /3 cos /3[(mi — 7712) cos^ 9\y + 1712 — ejm^Oj 



Afz["^2 — ^vy [/^ cos /3 + ejm^o sin /5] 



-tim^o [rh2 — ejm^o] [mi — eim^o] — M| cos^ — "^2) cos^ + m2 — ejm^o] 



(A.41^ 



M2[m2 — ejm^o] sin 6w[n cos /3 + ejm^o sin /3] 
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(A.42) 
(A.43) 



and the relation between Zj\f and is given by Eq. A. 24. 
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Appendix B 



In this appendix we list the imaginary parts which arise from the SUSY couphngs and were 
defined by the C^'s in Eqs. P^ - |55| , Eqs. and Eqs. {k stands for the corresponding diagram 

in Fig. 2, i.e. k = a,b,c,d,a',b',c',d',e, f,g,h). These contain the necessary CP-odd SUSY phases 
needed for both the production cross-section asymmetry, Ap, and the decay PRA in the top decay, 
Ad. 



ol ^ \zmHm{c} , (B.i; 

Ol ^ \Zl^\Hm{z^ZL-K+} , (B.4 

Ot = \Zl^\Hm{z^*L-K~} , {B.7 

0\ = -\m[K-M^] , (B.8 

Ol = lm{K-M^} , (B.9 

Ol = Im{ir+M^} , (B.IO 

O^ = -Im{ir+M^} , (B.ll 

Ol ^ -\\Zlf\Zl,"\'lm{L+L-*} , (B.I2: 

Ol, ^ \Zl''\Hm{z-^L~*K^*} , (B.13: 

Ol, ^ \Z}"\Hm{z^^L-*K-*} , (B.14 

Ol, ^ \Zlflra{zt*L^K-*} , (B.IS; 

Ol ^ \Zlflm{zt*L^K+*} , (B.ie: 

Ol ^ i=|Z^flm{z+:L+M2} , (B.I?: 

Ol ^ -^\Zlflm{zt:,L^M'} , (B.IS: 
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O] = -^Z\^\\v^{z^^L-*M^} , (B.20) 
^ -^|Z]°|2lm {Zr^L-*M2} , (B.21) 
Ol ^ -^\ZlW\'lm{L-L^Z^*Zt*} , (B.22) 

In order to avoid any confusions we chose the following notation for the indexes of different squark 
mass eigenstates: l,o = 1,2 are the indices for the two superpartners of the t, b, u, d, respectively, 
m = 1, 2 is the index for the two mass eigenstates of the charginos and n = 1 — 4 is the index for the 
four neutralinos mass eigenstates. Also L'^^K^ and K~ are defined in Eqs. A.9-A.11 and: 



S = Zf'Zf , (B.24) 



V^Miysin/? 
^ 1 I ( rrit 



sf2 \ VMvKsin/3 



j \Zf\^Zi^Zt^ -^\Z]fZt^L'^^ , (B.25) 



- (7=2 ■ ^^-''^ 
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Appendix C 



In this appendix we give the coefficients , D'l and Cq , Dq which appear in Eqs. |22|-^ Eqs. 



51| and Eqs. ^5]-^ {k stands for the corresponding diagram in Fig. 2, i.e. k = a,b, c, d, a', b', c', d', e, /, g 



and h and p=l,2,q=l — 7). These coefficients which are functions of masses and momenta are 
defined by the one- loop momentum integrals as follows p2|: 



d'^k 



1- h ■ h h 



m"^ [k^ — rrii] [{k + pi^ — 1712] [{k 



PzY - ml] 



(B.l) 



Do;D^;D^^(mi,m2,m3,m4,Pi,p2,P3>P4, (P1+P2) ,(^2+^3) ) 



d'^k 



1- h ■ h h 



in^ [fc^ — mf] [(k + pi)"^ — 1712] [{k + pi + P2Y ~ ^3] [(^ ~ ^4)^ ~ '^l] 



21 ' 



(B.2) 



where J2iPi = ^ is to be understood above. 



The coefficients are then defined through the following relations |33]: 



c 

D 



= Plf,PluC21 + P2MP2i.C22 + {PlP2}^„.C23 + g,,v^2i , 
■■ J5l/,Dii + P2mDi2 + P3mDi3 , 

= PIMPI1/D21 + P2piP2uQ22 + P3AtP3i/D23 + {P1P2} ^u^2A 
+ {Pm}iiv^2^ + {P2P?,} i,v^2& + fi'/.i.D27 , 



(B.3) 
(B.4) 
(B.5) 

(B.6) 



where {ab}^^ = a^b^, + a^b^. 

With the above definitions and notation, the coefficients Cq 's and C^^'s for the triangle diagrams 
in the production amplitude which appear in Eqs. ^2] - ^ , are given by: 



Im |Co; Cpqi^/ztfj, 



^O'l ^pq = (Co; Cpq 



)} 



(B.7) 
(B.8) 
(B.9) 
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Co^; C'^^ = C;^im-,^ ^ m,^,m^o ^ m^J , (B.IO) 

C°';C;; = Im{Co;Cpq(m^,m|,m|_^,m2,s,m^)} , (B.ll) 

Cj; = Co"'; C;;(mG - m^o ) , (B.12) 

'^o'sC'pg = Im{Co;Cpq(m|^,m|^,m|o,m2,s,m2)} , (B.13) 



The -Dq's ^'^d -Dpq's, for the box diagrams in the production amphtudes which appear in Eqs. ^^pT 
are given by: 



^o;^pg = Im{Do;Dpq(m|o,m|,m|^,m|,m2,m^,m^,m2,S,t)} , (B.15) 

^0^; = Dl- Dllm,^ m,^^m,, ^ ^ J , (B.16) 

^o; ^pg = ^o; Dpgiruu, rn^^,mu ^ ma, i ^ u) , (B.17) 

^0 ; = Dl- Dlgimi^ mr , mj^ m^, , m^o ^ m^^^) , (B.18) 



and for the top decay amphtudes, the coefficients Cg's and C^^'s which appear in Eqs. p5| - |6^ , are 
given by: 



Co ; C^g = Im {Co; Cpq(m|, , m| , m^, m^, m^^ ml) } , 



Cq; Cp„ — Cq ; Cp {rriG 



Im { Co ; Cpq (m|o , , m| , , , ) } 
Co; Cpg = Co'; C;g(mi , m^-^, m^o ^ m^„) . 



(B.19) 
(B.20) 
(B.21) 
(B.22) 



For the numerical evaluation of the above form factors see [^. 
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Figure Captions 



Fig. 1: The tree- level Feynman diagram contributing to ud — > tb. 

Fig. 2: The CP-violating, SUSY induced, one-loop diagrams for the processes ud — > tb and t — > 
W+b. 

Fig. 3: The allowed regions in the sina^ — sina^ plane for the NEDM not to exceed a) 1.1 x 10^^^ 
e-cm and b) 3 x 10^^^ e-cm. = 400 GeV and mo = 500 GeV is used. The shaded areas 
indicate the allowed regions. 

Fig. 4: The SUSY induced partially integrated cross-section asymmetry as a function of for 
Ms = 400 GeV, mi = 50 GeV and for = 2 TeV. With a) tan/3 = 1.5 and b) tan/3 = 35. 
See Ref. 24. 

Fig. 5: The SUSY induced partially integrated cross-section asymmetry as a function of mi, for 
several values of /i, M, = 400 GeV, mc = 450 GeV and for = 2 TeV. With a) tan/3 = 1.5 
and b) tan /3 = 35. 

Fig. 6: The SUSY induced partially integrated cross-section asymmetry as a function of mo, for 
several values of Ms = 400 GeV, mi = 50 GeV and for = 2 TeV. With a) tan /3 = 1.5 
and b) tan P — 35. 

Fig. 7: The SUSY induced partially integrated cross-section asymmetry as a function of tan (3, for 
several values of /x, M, = 400 GeV, mi = 50 GeV and for = 2 TeV. With a) mc = 350 
GeV and b) ma = 350 GeV. 

Fig. 8: The SUSY induced partial rate asymmetry in the top decays as a function of = 400 

GeV, ma = 350 GeV and for ^ = 2 TeV. With a) tan /3 = 1.5 and b) tan /3 = 35. See Ref. 24. 

Fig. 9: The SUSY induced partial rate asymmetry in the top decays as a function of mi, for several 
values of fi, Ms = 400 GeV and for = 2 TeV. With a) tan /3 = 1.5 and n = -160 GeV, and 
b) tan/3 = 35 and = 230 GeV. 

Fig. 10: The SUSY induced partial rate asymmetry in the top decays as a function of mc, for 
several values of /x, Ms = 400 GeV, mi = 50 GeV and for = 2 TeV. With a) tan/3 = 1.5 
and b) tan /3 = 35. 
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11: The SUSY induced partial rate asymmetry in the top decays as a function of tan (3, for 
several values of /x, M, = 400 GeV, mi = 50 GeV and for ^ = 2 TeV. With a) ma = 350 
GeV and b) ma = 350 GeV. 
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Table 1: The charginos and neutralinos masses for different values of {tan/3, mc}- The entries in 
the table are rounded to the nearest integer. 
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